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Abstract 

- > 


Scattering  of  elastic  waves  by  surfacetbreaking  planar  and  nontplanar  (branched)  cracks  ^ 

.•  /  As 


has  been  studied  in  this  paper.  Attention  has  been  focused  on  the  near?fleld  surface  displace- 
ments  and  the  crackjtip(s)  stressiintensity  factors.  For  planar  normal  cracks  the  stress! intensity 
factors  are  shown  to  agree  with  earlier  results.  Numerical  results  showing  normalized  vertical 
surface  displacements  are  presented  for  incident  body  and  surface  waves.  It  is  shown  that  the 
results  for  planar  and  branched  cracks  can  be  significantly  different  in  some  instances. 


I 


T 


Introduction 


Problems  of  elastic  wave  scattering  by  surface-breaking  and  near-surface  cracks  are  ^f  con¬ 
siderable  current  interest  for  ultrasonic  nondestructive  evaluatio^' Ultrasonic  scattering  by 
planar  cracks  near  or  at  the  free  surface  of  a  semi-infinite  elastic  homogeneous  medium  has  been 
studied  theoretically  by  many  authors.  References  to  recent  papers  on  this  subject  can  be  found 
in  1  ,  21  and  3  .  Some  experimental  works  on  surface-breaking  normal  planar  cracks  have  also 
appeared  (  4  -  6  ). 


Ultrasonic  scattering  by  surface- breaking  planar  and  branched  cracks  of  arbitrary  orienta¬ 
tion  is  the  subject  of  this  investigation.  To  our  knowledge  this  problem  has  not  received  much 
attention  in  the  literature.  An  approximate  solution  that  is  valid  at  low  frequencies  was 
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presented  in  [7]  for  SH  wave  diffraction  by  a  canted  surface-breaking  planar  crack.  Subse¬ 
quently,  a  hybrid  finite  element  and  eigenfunction  technique  was  used  in  [8]  to  study  SH  wave 
diffraction  by  planar  surface-breaking  canted  crack. 

In  this  paper  we  use  the  same  hybrid  technique  as  in  [8]  to  study  the  scattering  of  in-plane 
body  and  surface  waves  by  canted  planar  and  normal  surface- breaking  branched  cracks.  We 
focus  our  attention  to  the  near-field.  Numerical  results  are  presented  for  the  vertical  surface 
displacement  amplitudes  near  the  base  of  the  crack  and  crack-tip(s)  stress-intensity  factors. 

Formulation  and  Solution 

Consider  a  homogeneous,  isotropic,  and  linearly  elastic  medium  with  a  surface-breaking 
crack  of  arbitrary  orientation  and  shape  as  shown  in  Figure  1.  Assume  the  displacement 
u(*  ,y  ,t )  at  a  point  P  to  be  time-harmonic  of  the  form  u(x,y)e~ twt  where  is  the  circular 
frequency.  Then  u  satisfies  the  equation  of  motion  in  y  >  0  (at  points  not  on  the  crack) 

+  (A  +  n)  vv- «  =  -  (1) 

where  A,  p  are  Lame  constants,  p  the  mass  density  and  the  factor  e  -,w<  has  been  dropped. 

The  solution  of  (1)  can  be  expressed  in  terms  of  longitudinal  and  shear  wave  potentials,  d> 


and  0,  in  the  form 


«  -  v*1  +  v  *  ) 


(2) 


F urthermore,  in  a  homogeneous  half-space,  ^  and  tl>  can  be  expressed  in  infinite  series  of  mul¬ 
tipolar  potentials  as  {9j. 


*-  E  h  +S  +  KK 


*=  E  (an^nP  +  6n<| 


(3) 
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where  expressions  for  ip£,  4>  and  ip*,  can  be  found  in  [9].  The  coefficients  an ,  bn  are 
found  by  satisfying  the  appropriate  boundary  conditions. 

The  representation  (3)  is  not  useful  for  satisfying  the  boundary  conditions  on  the  crack 
surface.  For  this  reason,  a  different  representation  is  needed  in  this  near-field  region.  In  this 
paper,  the  region  inside  the  fictitious  boundary  B  (Fig.  1)  is  divided  into  finite  elements  having 
Nj  number  of  interior  nodes  and  Ng  number  of  boundary  nodes. 

For  the  finite  element  representation  in  region  II,  the  energy  functional  is  taken  to  be 


F  *  —  J  J*  (o-e*  -  pw2u*u*J  dzdy 

*n 

~  \  £  [?B  '-B  +  Fb'Hb)** 


(4) 


where  indicates  complex  conjugate  and  o  and  e  are  column  vectors  defined  as 


!m  <f)m  (<«•%„’**„) 


Superscript  ‘T’  denotes  transpose.  The  ?B  and  Ug  represent  the  traction  and  displacement  on 
B ,  respectively. 

It  is  assumed  that  the  displacement  field  within  the  j1*1  element  is  represented  in  terms  of 
the  shape  functions  L-{x,  y)  and  elemental  nodal  displacements  {?  ?}  as 

=  £  Ljt} 

y»i 

where  each  q  j  has  two  components  and  along  the  x  and  y  directions,  respectively. 
The  .V  represents  the  number  of  nodes  in  each  element. 
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The  o  *•  and  ef*-  are  computed  by  substituting  Eq.  (7)  into  strain-displacement  relations 
and  these,  in  turn,  into  the  stress-strain  relations.  Using  these  in  Eq.  (4),  we  get 


F  *  qfT  Sjj  qj  +  q*T  SIB  qg  +  qgT  S Bih 


(8) 


„  „*T  p  (1)  _  p*T(l)  „ 

+  lB  bBB  iB  ~  i  -B  IB 

in  which  qj  =  qj^\  qg  -  q^\  P^  =  P^  and  the  elemental  impedance  matrices  5f- ■  are 


defined  as 


[S<]=  I  {  {[B'fmB*]-  Py\L\T\L\)  dzdy  . 


(9) 


In  Eq.  (9), 


l  B*)~ 


o 

0 

dx 

n 

_a_ 

U 

By 

0 

_d_ 

By 

dx 

L  y  0  L  2 
0  Lx  0 


WI 


Note  that  [LJ  is  a  2  x  2 Ng  matrix. 

For  an  isotropic  material  [D]  is  given  by 


Ae  +  2/i  e  Ag  0 

A-  Ag  +  2ne  0 

0  K 


e 

0 


where  Ag  and  P  3  are  the  Lame’s  constant. 

To  find  the  constants  an  ,  6n  appearing  in  (3)  and  the  nodal  displacements  in  region  II,  it 
is  necessary  to  use  the  continuity  of  displacement  and  traction  on  B.  This  is  discussed  in  the 
following. 
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The  incident  displacement  fields  will  be  assumed  to  arise  from  the  incident  plane  P  and  SV 
waves,  and  their  reflections  from  the  free  surface  y  =  0.  The  case  of  incident  Rayleigh  waves 
will  also  be  considered. 

Let  us  suppose  that  in  the  absence  of  the  crack  the  free  field  is  the  sum  of  the  incident  and 
reflected  fields,  that  is 

UJ0)  m  „*>  +  uJr )  (y  =  1>2)  .  (10) 

For  the  Rayleigh  wave  u  is  the  associated  displacement. 

The  total  field  outside  B  then  is 

«/  ■■/*)  +  U  *  1.2)  (11) 

where  l  is  given  by  (2)  and  (3). 

Using  (2)  and  (3),  the  displacements  at  the  nodes  on  B  can  be  written  as 

*  lG){?>  (12) 

where  (G  ]  is  a  2Ng  x  2Ng  matrix  formulated  in  Appendix  A  and  vector  {a}  is 

{?}* 

Similarly,  using  (2)  and  (3)  in  the  stress  strain  relation,  the  traction  at  the  nodes  on  B  can 
be  expressed  in  the  form 

(13) 

where  [Fj  is  also  a  IN g  x  2 Ng  matrix  defined  in  Appendix  A. 

To  express  (o^1 )}  in  terms  of  {qjf  )},  we  use  the  expression  for  the  virtual  work  done  on 
the  boundary  B  ,  which  is 

**  -  f  {«?s(1)>T  (14) 

B 
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where  superscript  (1)  denotes  the  total  field  in  region  I  (outside  B). 

Because  of  the  continuity  of  displacements  and  traction  on  B ,  we  have 


tl"  -  d"  =  iL0)  +  il‘> 

(15) 

tP-tP-tP+tlr' 

(16) 

where  superscript  (2)  denotes  the  total  field  in  region  11. 

Substituting  (12),  (13),  (15),  and  (16)  in  Eqn.  (14),  and 

noting  that  Sq ^  =  6q^\  we 

obtain  from  Eqn.  (14) 

Sr  =  {Sa*}T{P^} 

(17) 

where  /*£*)  is  given  by 

«#’>  -  [*!<«> + {?r> 

(18) 

Here 

\n\-  ^\c*\T\F]dr 

(19) 

and 

(20) 

Equations  (19)  and  (20)  are  approximated  by 

[It]  =  {G*]T[F]*RM 

and 

(21) 

{P^}  =  G*\T 

(22) 

where  R  A0  is  the  arc  length  between  two  adjacent  boundary  nodes  on  contour  B .  Note  that 

the  first  two  rows  and  last  two  rows  of  Eqs.  (21)  and  (22)  are 

multiplied  by  R  instead  of 

R  &0,  because  they  correspond  to  the  first  and  last  boundary  nodes,  respectively. 

Substituting  Eqn.  (12)  in  Eqn.  (8)  and  taking  the  variation,  we  obtain  a  set  of  simultane¬ 
ous  equations  which  may  be  written  in  matrix  form  as 


S!I 

S1B  GI 

il 

/>  *  T  c  T 
G  SIB 

g*tsbbg 

a 

-G*TsBBili0)  +?il). 


(23) 


The  first  equation  of  Eqn.  (23)  can  be  written  as: 

?/  “  ~  SIJ  1  [5/B  +  SIB  fir0^]  ■  (24) 

The  second  equation  can  be  written  as 

G*TS$,U+G*TSbbG  *--G*TSBBij?)  +#)  .  (25) 

Substituting  Eqs.  (18)  and  (24)  into  Eq.  (25),  we  obtain 

[G  *  T  “  SIB  SlFl  SIB  j  G  *  ~  *]  W  * 

(26) 

-  G*T  [sgg  - 


In  Eq.  (26),  the  generalized  coordinates  (a)  are  the  only  unknowns.  Therefore,  {a}  can  be 
evaluated.  Once  (a)  are  known,  the  near  and  far  displacement  and  stress  fields  can  be  deter¬ 
mined. 


Numerical  Results  and  Discussion 

In  this  paper,  the  boundary  B  enclosing  the  interior  region  is  not  a  complete  circle,  and  so 
the  potentials  <t>  and  <j>*,  tf>‘  cannot  be  expanded  in  cylindrical  wave  functions  as  was 

done  in  [1],  So  the  integrals  giving  these  potentials  and  their  derivatives  were  evaluated  numeri¬ 
cally  for  every  node  on  B  .  The  details  are  discussed  in  (10}. 


The  hybrid  method  is  employed  to  study  scattering  by  P,  SV,  and  Rayleigh  waves  by 


three  types  of  surface  breaking  cracks:  a  vertical  crack  (Fig.  2  with  a  =  90  ’ ),  a'  45 "  inclined 
crack  (Fig.  2,  with  a  =  45"),  and  a  vertical  branched  (Y)  crack  (Fig.  3).  The  Poisson  s  ratio  of 
the  material  is  taken  to  be  1/3. 

Stress  intensity  factors  at  the  tips  of  the  cracks  were  calculated  and  for  the  particular  case 
of  a  planar  surface-breaking  normal  crack  they  were  found  to  agree  well  with  the  results  of  [1]. 
These  are  shown  in  Figure  4. 

Next,  normalized  stress-intensity  factors  are  shown  in  Figures  5-7  for  incident  P,  SV,  and 
Rayleigh  waves.  It  is  seen  that  for  P  and  Rayleigh  waves  the  stress-intensity  factors  at  the 
crack-tips  of  the  three  types  of  cracks  are  quite  different,  particularly  at  high  frequencies.  This 
is  particularly  significant  for  the  branched  crack,  even  though  the  branches  are  quite  small. 

The  surface  displacements  at  y  =  0  are  calculated  by  using  (3)  in  (2)  after  {a  }  are  calcu¬ 
lated.  Normalized  values  of  u  )  are  presented  in  Figures  8-12.  For  each  type  of  crack  men¬ 
tioned  above  five  cases  of  incident  waves  were  considered:  plane  P  wave  incident  at  0  *  and  45  ° , 
plane  SV  wave  incident  at  0  *  and  45  * ,  and  Finally  Rayleigh  wave.  Some  representative  results 
are  shown  here. 

Figures  8  and  9  show  the  scattered  vertical  surface  displacement  amplitudes  for  a  Rayleigh 
wave  incident  from  the  left  on  a  normal  planar  and  branched  crack.  It  is  seen  that  there  are 
large  differences  in  the  forward  direction  between  the  two  cases  as  the  frequency  becomes  large. 
In  the  backward  direction,  however,  the  differences  are  not  very  significant.  Figures  10  and  11 
show  the  results  for  an  incident  SV  wave  moving  vertically  as  well  as  at  45  *  to  the  vertical. 
Large  differences  are  found  for  vertical  incidence,  but  not  in  the  other  case.  Finally,  in  Figure 
12  is  shown  the  case  of  a  Rayleigh  wave  incident  on  a  canted  crack.  This  figure  is  to  be  con¬ 
trasted  with  Fig.  8.  The  large  contrast  shown  clearly  distinguishes  a  canted  crack  from  a  nor¬ 


mal  crack. 
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Concluaion 

Model  calculations  of  elastic  wave  scattering  by  surface-breaking  planar  and  non-planar 
cracks  have  been  presented.  These  calculations  show  that  near-field  surface  displacements  due 
to  scattering  by  planar  and  branched  cracks  are  quite  different  even  when  the  branches  are 
small.  Also,  it  is  found  that  signatures  of  normal  and  canted  cracks  are  very  dissimilar.  These 
characteristic  differences  can  be  used  to  discriminate  between  the  various  cases. 

Although  the  results  presented  here  are  for  homogeneous  medium,  the  technique  can  be 
generalized  to  study  cracks  in  a  composite  medium.  These  are  presently  under  investigation 
and  will  be  reported  later. 
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APPENDIX  A 

Formulation  of  Matrices  [G]  and  [F] 

As  mentioned  before,  the  scattered  nodal  displacement  vector,  {qjfh i  was  formed  by 
evaluating  )  and  e  )  at  Ng  number  of  points  on  contour  B .  Thus,  we  have: 

{«£>)-  (C)  {;}  (A-1) 

where, 


}> = Kbj’  •  •  •  *  \bn  •  v,’  ■  ■  ■  •  %bn  v 


(A-2) 


{?}={«!>  *  •  •  .  »  *  l>  ‘  ‘  •  .  bN„ 


(A-3) 


If  [C]  is  partitioned  as, 


GXA  j  GXB 

I 

GYA  |  GYB 


2  Ng  X  2Ng 


[A-4] 


then  each  of  the  x  Nq  submatrix  can  be  evaluated  from  Eqns.  (2)  and  (3)  at  )  on  B 


as 


-  «,+<,) 

(GAB),,,  -  (*/,,  +<,) 

(GYA  ),„  -  «„-<,) 


(A-5) 


(GYB),n  •  [*l,  -*.*,)  .  •  " 
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The  parameter  n  in  the  summation  series  of  Eqn.  (3)  is  taken  from  (- Ng  /2  —  l)  to  Ng  /2 


■'* 

A 

I 


I'i 

ll|! 

■S 

•it 


for  numerical  purposes.  Hence,  n  in  the  first  column  to  the  last  column  of  each  submatrix 
corresponds  to  {-Ng  2  -  1)  to  Ng/ 2,  respectively. 

To  formulate  the  matrix  F  ,  the  components  Tx  and  T  of  the  traction  vector  T  were 
calculated  at  each  nodal  point  on  B .  If  the  radius  vector  of  that  point  makes  an  angle  /?  with 
the  x-axes.  then 


J#  I 

i  Si 


M 

I 


8 

Si 


:*«» 

i# 

i  v 

.#rJ 


7*r  =  a  Ix  cos#  +  a  sin# 


•  -  *rvcos*  +  %ysi“*  • 


Evaluating  tr  c  J,  and  a  Jj  J  at  number  of  points  on  contour  B  and  substitut¬ 
ing  in  Eqn.  (A-6),  we  have  the  scattered  nodal  stress  vector,  {ajf  )},  as  in  Eqn.  (13), 

(f|«  (A-T) 


where 


<»£’>-  {t,b .  r. 


and  {a}  is  defined  in  Eqn.  (A-3) 
[F]  is  partitioned  as 


•M 


FA/4  |  FYAB 


FV04  I  FYB 


2AB  x  2JVB 


1 

1 


Each  of  the  ATg  x  /Vg  submatrix  can  be  evaluated  at  (zf-,  yf- )  on  B  as. 
(FA‘‘4  ).n  =  [(A  +  2^n,zx  +  + n,x„ )  +  M^yy  "  * rf.zy)]  cos^ 


+  u(2<t>  P  ^  P  )  sin/? 

n,ry  ,yy  rn,xx'  u 


■*_» 

$ 


CO'  .1 


»V* 

vV 


r-ft 

voM 

V  f 


m 
% 

YM 


;$| 

*.»>M 

■  ,V, 


Vo 


>  V» 

w 

'.W 

V'I|- * 

W:t. 


«•$ 

$4 

$ 


.*,v] 


(™)m  =  [(A  -2^(0/ 


•1  yj  -n.ry 


—  u( 2o  ^  —  v  ~  -  v  „  )  sinJ 

^  n  ,zy  n  .yy  m 


(FK4  )•  =  m(2o  P  -  v  } 

v  •  tn  n  ,ry  n 


xy  n  .yy  #»  x: 


cos  J  - 


[(A  —  2fi){6  P -  vf  zy)  -  MonPIZ  -  enPtJ,)]  sin/? 

(™>«.  =^n,ty  *  -«•/.„)  coSJ  - 

[(A  ^  2M)«5w  -  ^nSiy)  +  *(*/.„  *  ^n5iy)]  sin/? 
Parameter  n  ranges  from  ~(Ng  /2  -  1)  to  jVg  /2  as  discused  before,  and  »  *  1  to  Ng 
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List  of  Fifuros 

1.  Geometry  of  s  mar- surface  inhomogeneity 

2.  Surface-breaking  canted  crack* 

3.  Surface-breaking  normal  branched  crack 

4.  Comparison  of  s  tress-intensity  factors  calculated  by  the  present  method  and  those  given  in 
[11].  Incident  wave  is  a  Rayleigh  wave  propagating  in  the  x-direction. 

5.  Stress-intensity  factors  for  various  crack  geometries  due  to  incident  P  wave 

6.  Stress-intensity  factor  for  various  crack  geometries  due  to  incident  SV  wave 

7.  Stress-intensity  factor  for  various  crack  geometries  due  to  incident  Rayleigh  wave  pro¬ 
pagating  in  the  x-direction 

8.  Normalized  scattered  vertical  surface  displacement  amplitude  due  to  Rayleigh  wave 

incident  on  a  normal  surface- breaking  planar  crack 

9.  Normalized  scattered  vertical  surface  displacement  amplitude  due  to  Rayleigh  wave 

incident  on  a  normal  surface-breaking  branched  crack 

10.  Normalized  scattered  vertical  surface  displacement  amplitude  due  to  SV  wave  incident  on 
a  normal  surface-breaking  planar  crack 

11.  Normalized  scattered  vertical  surface  displacement  amplitude  due  to  SV  wave  incident  on 
a  normal  surface-breaking  branched  crack 

12.  Normalized  scattered  vertical  surface  displacement  amplitude  due  to  Rayleigh  wave 

incident  on  a  normal  surface-breaking  canted  crack 
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45°  INC.  CRACK 


